Probability of anomalously large Bit-Error-Rate in long haul optical transmission 
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We consider a linear model of optical pulse transmission through fiber with birefringent disorder 
in the presence of amplifier noise. Both disorder and noise are assumed to be weak, i.e. the average 
bit-error rate (BER) is small. The probability of rare violent events leading to the values of BER 
much larger than its typical value is estimated. We show that the probability distribution has a 
long algebraic tail. 
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Optical fibers are widely used for transmission of infor- 
mation. In an ideal case, information carried by pulses 
would be transmitted non-damaged. In reality, however, 
various impairments lead to the information loss. Noise 
generated by optical amplifiers and fiber birefringence are 
the two major impairments in high-speed fiber commu- 
nications. The amplifier noise is short-correlated in time, 
while the birefringence varies significantly along the op- 
tical line and is practically frozen in time, since the char- 
acteristic temporal scale of such variations is long com- 
pared to the signal propagation time through the entire 
fiber line. Coexistence of two different sources of ran- 
domness characterized by two well-separated time scales 
is pretty common in statistical physics of disordered sys- 
tems. A classical example would be the glassy behav- 
ior driven by short-correlated thermal noise in a system 
with frozen structural disorder (see, e.g., Q for review). 
An extremely non-Gaussian statistics of the observables 
is an important feature of the disordered/glassy system. 
In this letter we demonstrate that strong deviation from 
Gaussianity is also typical for optical fiber telecommuni- 
cation systems. 

Birefringent disorder is caused by weak random ellip- 
ticity of the fiber cross section. Birefringence splits the 
pulse into two polarization components and it also leads 
to pulse broadening [3- 0- Dl] - This effect known as po- 
larization mode dispersion (PMDJ_have been extensively 
studied experimentally 0, 0, Q, E3 EH and theoreti- 
cally [ill flil . PMD is usually characterized by the 
so-called PMD vector that was found to obey Gaussian 
statistics [lllll^.ll3j . It was also shown, e.g. in ^ij, that 
first order PMD compensation corresponding to cancella- 
tion of the PMD vector on the carrier frequency is exper- 
imentally implementable. Higher-order generalizations 
of the PMD vector (introduced to resolve complex fre- 
quency dependence of the PMD phenomenon with higher 
accuracy) as well as a suggestion on how to compensate 
for PMD in the higher orders have been also discussed 
and implemented experimentally, e.g. in EB.]. This 
works were focused on separating effects caused by PMD 



from other potential impairments. Common wisdom hid- 
ing behind this strategy says that one should start with 
evaluating the two impairments separately and then es- 
timate the joint effect on the optical telecommunication 
system performance taking them on equal footing. In this 
letter we challenge this equal-footing approach. We show 
that the effects of temporal noise and structural disorder 
on the overall system performance may not be separated 
since the bit-error-rate (BER) strongly depends on a re- 
alization of birefringent disorder. Thus, the probability 
distribution function (PDF) of BER and especially its 
tail corresponding to large values of BER are the objects 
of prime interest and practical importance for describing 
the probability of the system outage. 

Our letter is organized as follows. We start with dis- 
cussing the dynamical equations for an optical pulse evo- 
lution in a birefringent fiber also influenced by the am- 
plifier noise. BER produced by the amplifier noise for a 
given realization of the birefringent disorder is analyzed 
first. Then we describe the PDF of BER, found by av- 
eraging over many realizations of disorder. Some general 
remarks conclude the letter. 

The envelope of the electromagnetic field propagating 
through optical fiber in the linear regime (i.e. at rela- 
tively low pulse intensity), which is subject to PMD dis- 
tortion and amplifier noise, satisfies the following equa- 
tion HHE3 

d z & - iA(z)& - m(z)d t & - id(z)d%& = $(z, t), (1) 

where z is the position along the fiber, t is the retarded 
time (measured in the reference frame moving with the 
optical signal), £ is the amplifier noise, and d is the chro- 
matic dispersion coefficient. The envelope ^ is a two- 
component complex field where the components stand 
for two polarization states of the optical signal. Birefrin- 
gent disorder is characterized by two random Hermitian 
2x2 traceless matrices A and rh measuring fiber bire- 
fringence in the first- and second-order and correspond- 
ing to the first two terms of the expansion in uj — ujq 
with uj and loq being the signal and carrier frequencies 
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respectively. The disorder is frozen at least on all the 
propagation related time scales, i.e. the two matrices can 
be considered to be ^-independent. The random matrix, 
A(z), can be excluded from the consideration by means 
of the transformation to the coordinate system rotating 
together with polarization state of the signal at the car- 
rier frequency: \P — ► V\P, £ — » V£ and rh — » VrhV" 1 . 
Here, the unitary matrix V(z) is the ordered exponen- 
tial, Texp[j J* dz'A(z')] defined as the solution of the 
equation, d z V — iAV, with V(0) = 1. Hereafter we 
use the notations rh, £ and \P for the renormalized ob- 
jects whereas the original objects will be no more referred 
to. With this notation the equation for & has the same 
form as Eq. (JTJ if we set A — 0. The solution of the 
renormalized equation can be partitioned into a sum of 
a homogeneous contribution, insensitive to the additive 
noise and the inhomogeneous one, <f> and tp, respectively: 

V = W(*)*b(t), <t> = [dz'WWW-^z'Mz'rf, (2) 
Jo 



W(z) — exp 



* / dz'd(z')d? 





Texp 



dz'rh(z')dt 



, (3) 



where tf'o describes the input pulse shape (at z = 0). 

We assume the optical system length Z is much larger 
than the distance between the next-nearest amplifier sta- 
tions. (The stations are set to compensate for losses 
in the pulse intensity.) Coarse-graining on the inter- 
amplifier scale allows treating amplification in the con- 
tinuous limit. Zero in average additive noise, £, is the 
amplification leftover. The amplifier noise has Gaussian 
statistics 0] and its correlation time (set by quantum ex- 
citation processes in amplifiers) is much shorter then the 
pulse width. Therefore, ^-statistics is fully determined 
by its pair correlation function 

(&(*i,ti)Q(*2, t 2 )> = DtfcfS (zi - z 2 ) 8{tx - ta), (4) 

where the coefficient characterizes the noise strength. 

Averaging over birefringent disorder is of different na- 
ture: statistics here is collected over different fibers or, 
alternatively, over different states of the same fiber col- 
lected over time. (Birefringence is known to vary on a 
time scale essentially exceeding the pulse propagation 
time.) The matrix rh can be expanded in the Pauli matri- 
ces, rh(z) = hj(z)aj, where hj is a real three-component 
field. The field is zero in average and short-correlated in 
z. It enters the observables described by Eqs. ijL'l.'ill in an 
integral form and, according to the central limit theorem 
(see, e.g., 0]), can be treated as a Gaussian random 
field described by the following pair correlation function 



(hi(zi)hj(z2)) = D m 5ij5(zi ~ z 2 ), 



(5) 



where D m characterizes the disorder strength. Note that 
even though the original rh entering Eq. Q is anisotropic 
the isotropy of hj [implied in Eq. (JjJJ] is restored as a 
result of the rh — > VrhV^ 1 transformation. 



We discuss here the case of the so-called return-to-zero 
(RZ) modulation format when pulses in a given frequency 
channel are well separated in t. Detection of a pulse 
at the fiber output corresponding to z = Z requires a 
measurement of the pulse intensity, /, 



1=1 dtG(t) \K¥(Z, t)\ 



(6) 



where the function G(t) is a convolution of the electrical 
(current) filter function with the sampling window func- 
tion (limiting the information slot). The linear operator 
KL in Eq. JfjJ) stands for an optical filter and may also in- 
corporate a variety of engineering "tricks" applied to the 
output signal \P(Z,i). Ideally, I takes a distinct value if 
the bit is "1" and is negligible if the bit is "0" . Both the 
noise and the disorder enforce I to deviate from its ideal 
value. One declares the output signal to code or 1 if 
the value of I is less or larger than the decision thresh- 
old 7o- The information is lost if the output value of the 
bit differs from the input one. The probability of such 
event should be small (this is a mandatory condition for a 
successful fiber line performance) i.e. both impairments 
typically cause only small distortion to a pulse. Formally, 
this means: D^Z <C 1, D m Z <C 1, where both the initial 
signal width and its amplitude are rescaled to unity. 

Below we focus our analysis on the initially "1" (uni- 
tary) bit. We do not consider the evolution of a zero 
bit since it does not contribute to the anomalously large 
values of BER, which we are mainly interested to de- 
scribe. The probability to loose the unitary bit is B = 
J /o dIP {I). Here P(I) is the PDF of the signal inten- 
sity © (that fluctuates due to the noise £) for the initial 
signal corresponding to the bit "1" . In engineering 
practice BER is measured collecting the statistics over 
many initially identical pulses. (As different pulses sense 
different realizations of the noise, this averaging over dif- 
ferent pulses is actually equivalent to the noise averag- 
ing.) Repeating the measurement of B many times (each 
separated from the previous one by a time interval larger 
than the characteristic time of the disorder variations or, 
alternatively, making measurements on different fibers) 
one constructs the PDF of B, S(B). The PDF achieves 
its maximum at Bq, a typical value of B. Even though 
the average distortion of a pulse caused by the noise and 
disorder is weak, rare but violent events may substan- 
tially affect the optical system performance. The prob- 
ability of such rare events is determined by S(B) taken 
at B 3> Bq. Our further analysis is focused on this PDF 
tail. 

Mentioning the importance of accounting for an ap- 
propriate form of the linear operator K., we restrict our 
discussion to two possibilities, when the optical filtering 
is accompanied by an overall time shift (this operation 
is usually called "setting the clock") or by a compen- 
sation achieved by insertion of an additional piece of 
fiber with adjustable birefringence. Optical filter is re- 
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quired to separate different frequency channels. In addi- 
tion the filter smoothes out an otherwise strong impact 
on the pulse caused by amplifier noise temporal ultra- 
locality. "Setting the clock" procedure is formalized as, 
K, c i^ — \P{t — t c {), where t c i is an optimal time delay. 
The so-called first-order compensation means 



= exp 



dz rh{z)dt 



(7) 



where the usual exponent (instead of the ordered one) 
enters Eq. ©. In the engineering practice or laboratory 
experiments all the three strategies (along with some ad- 
ditional others) can be applied simultaneously. Below we 
imply that the optical filter is always inserted. As far as 
the "setting the clock" and the first-order compensation 
procedures are concerned we intend to compare those ef- 
fects with the pure (no compensation) case. 

The first step in our calculations is to find the value 
of B for a given realization of disorder. In this case 
the yj-part of the envelope & is a constant (depen- 
dent on the disorder) while the ^-contribution fluctu- 
ates. One obtains from Eqs. \\'2\ \\ that <p is a zero 
mean Gaussian variable, with the pair correlation func- 
tion {4> a {z, ti)<f)*p{z, £2)) = D^z5 a p5(ti — is). Note that 
statistics of <j> is sensitive neither to the chromatic disper- 
sion coefficient d no to the birefringence matrix to. Thus, 
averaging over the noise statistics is reduced to a Gaus- 
sian path integral over cj). The inequality D^Z <§C 1 jus- 
tifies the saddle-point evaluation of the integral and also 
allows to estimate \nB as lnP(Jo). Thus, one finds that 
the product D^ZhiB is a negative quantity of order one, 
insensitive to the noise characteristics. This quantity de- 
pends on the initial signal profile disorder hj, integral 

chromatic dispersion J Q Z d{z)dz, and also on the details 
of the measurement and compensation procedures. 

We next analyze the dependence of B on the birefrin- 
gent disorder. The smallness of D^Z means that even 
weak disorder can produce large deviations in B from 
its typical value B , \nB ~ — (D^Z)^ 1 , corresponding 
to hj = 0. Therefore \n(B/Bo) can be represented as 
a series in hj. The leading contribution is found upon 
expanding the ordered exponential © [entering ip in ac- 
cordance with Eq. in a series in h, followed by sub- 
stituting the result for cp into Eq. |JB}, and perform- 
ing the saddle-point calculations, aiming to find P{Io)- 
The smallness of D^Z enables us to substitute lni? by 
lnP(io). Then, in the second-order in h, one obtains 

HB/B ) = T/(D t Z), T = Ml tf 3 + (8) 

+/4 / dzi / dz 2 [hi(zi)h 2 (z 2 )-h 2 {z 1 )hi(z 2 )], 
Jo Jo 

where Hj = J dz hj(z) and the initial pulse l^b is as- 
sumed to be linearly polarized. The coefficients /j, enter- 
ing Eq. © are sensitive to the initial signal profile 



chromatic dispersion and to the measurement procedure. 
The coefficients can be calculated only numerically, how- 
ever, some conclusions based on the general form of Eq. 
(jHJ) can be drown without this detailed calculations. 

If no correction procedure is applied to the output 
signal then \ii ~ 1, whereas [i\ and fi' 2 are corrections 
related to the temporal asymmetry of the pulse (pro- 
duced by the optical filter) and to the so-called the 
pulse chirp (generated by the dispersion term), respec- 
tively. If no compensation is applied the leading role is 
played by the first term in Eq. JSJ), and one arrives at 
InP w -Y' 2 j{2[i\D m Z). If the "setting the clock" pro- 
cedure applies (with t c i dependent on h: t c i = to — H3, 
with to being the optimal time shift without disorder), 
the first-order term in the right-hand side of Eq. (jHJ 
cancels out and the second-order term is dominated by 
^{Hl + H" 2 ). According to Eqs. I|2I3I) the statistics of B 
depends on the integral dispersion J Q Z dzd(z). In practice 
some special technical efforts (dispersion compensation) 
are made to reduce it, so that usually this quantity is 
small. If the initial pulse is real (there is no chirp) and if 
the integral dispersion is small, then n' 2 is negligible. In 
this case Eq. (JjJJ leads to the following PDF of T 



V{T) = 



1 



2fi 2 D m Z 



exp 



2fi 2 D m Z 



(9) 



If the first-order compensation scheme described by Eq. 
is used the first two terms in the right-hand side of 
Eq. JHJ cancel out, /ii = [i 2 = 0, and if in addition [i' 2 ^ 
one arrives at 



P(T) = (2D m Z^ 2 )- 1 cosh- 1 ( 



71T 



(10) 



replacing Eq. Comparing Eqs. (|9I10() with the re- 

sult for the no-compensation case we conclude that the 
compensations lead to much narrower PDF. If fj,' 2 is also 
(i.e. if the first order compensation is applied and the 
output signal is not chirped) then higher-order terms in 
the expansion of T in hj should be accounted for. This 
case will be discussed elsewhere. 

Finally, we are in a position to analyze the PDF of 
B, S(B). Expressing T via B according to © followed 
by substituting the result into Eq. © for the PDF of T 
leads to the following expression for the remote (B ^> Bq) 
tail of the PDF of B 



S(B)dB 



B$ dB 



D, 



2fi 2 D„ 



(11) 



The applicability range for Eq. Hll|) is given by 1 ^ 
r > D m Z that transforms into 1/(D^Z) > ln(B/B ) > 
D m /D^. If the first-order compensation is applied Eq. 
© should be replaced by Eq. Hl()(l . thus leading to the 
same expression (jl 1|> for the tail with /j, 2 replaced by 
H 2 /tt. The dependence of the PDF on BER is illustrated 
in Fig. n One also finds from Eq. that the so-called 
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FIG. 1: Schematic log-log plot of the PDF of Bit-Error-Rate. 

outage probability, O = J B dBS(B), where f?* is some 
fixed value taken to be much larger than Bq, is estimated 

as InO~(Uf/U m )ln(flo/B*). 

There is also a universal remote tail of S(B) corre- 
sponding to huge fluctuations of the disorder when the 
signal is almost destroyed by the fluctuations. In this 
extreme case / is close to the threshold value Iq already 
at £ = 0. Then B is nothing but a probability of such an 
event (i.e. the noise configuration £) that the inhomoge- 
neous contribution <p would fill the remaining small gap 
between / and Iq. Thus one gets \nB ~ —cj) 2 /(D^Z). 
The value of 4> is proportional to the deviation Sh of the 
/i-disorder term from such a special configuration, h sp , 
which gives I = Iq at £ = 0. Therefore, one estimates, 
ln£> ~ — (Sh) 2 Z/D^. The logarithm of the probability for 
this disorder configuration to occur is a sum of two con- 
tributions: ~ —l/(D m Z), corresponding to the special 
configuration, h sp , and the other one, ~ 5h/D m , corre- 
sponding to Sh. One arrives at the following expression 
for the probability density 

where C\ and C2 are constants of order one. Eq. Q12|l 
holds when (D m Z) 2 < D i Z\hxB\ < 1. Note, that the 
remote tail 11120 decays with B faster than the algebraic 
one l|TT|). 

Eq. iffi"]) describes our major result: The PDF of BER 
has a long algebraic tail. The exponent a of the alge- 
braic decay is proportional to the ratio of the amplifier 
noise variance, D^, to the birefringent disorder variance, 
D m . This statement clearly shows that effects of noise 
and disorder are, actually, inseparable. Another inter- 
esting feature of Eq. (jl 1|> is that the exponent a is Z- 
independent. The only Z-dependent factor in the final 
result |(TJ| is the overall normalization factor Bq . Note 
also that even though an extensive experimental (labo- 
ratory and/or field trial) confirmation of this statement 
would be very welcome, some numerical results, consis- 



tent with Eq. I|ll|) are already available. Thus, Fig. 2a 
of [23 re-plotted in log-log variables shows a linear rela- 
tion between In S and InB, i.e. just the algebraic decay 
predicted by Eq. l|TTjl . 

Summarizing, this letter is the first brief report on the 
new method/approach. We plan to discuss many related 
but more specific issues, e.g. more complex compensation 
strategies (periodic and higher-order compensation) in 
subsequent publications in specialized optical journals. 
However, it is important to stress that the major result 
of the paper that the PDF of BER is characterized by a 
long tail (algebraic or faster, but slower than exponential) 
is universal. 
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